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Frustration, or the competition between interacting components of a network, is often responsible for tlie 
complexity of many body systems, from social fll and neural |2| networks to protein folding |3| and magnetism 
l4H6l. In quantum magnetic systems, frustration arises naturally from competing spin-spin interactions given by 
the geometry of the spin lattice or by the presence of long-range antiferromagnetic couplings. Frustrated mag- 
netism is a hallmark of poorly understood systems such as quantum spin liquids, spin glasses 1 7 , 8 1 and spin ices 
|91 , whose ground states are massively degenerate and can carry high degrees of quantum entanglement IIOIII II . 
The controlled study of frustrated magnetism in materials is hampered by short dynamical time scales and the 
presence of impurities, while numerical modeling is generally intractable when dealing with dynamics beyond 
A*' ~ 30 particles 1 12 |. Alternatively, a quantum simulator 1 13 1 can be exploited to directly engineer prescribed 
frustrated interactions between controlled quantum systems, and several small-scale experiments have moved 
in this direction 1 11 14 - 19|. In this article, we perform a quantum simulation of a long-range antiferromag- 
netic quantum Ising model with a transverse field, on a crystal of up to A'^ — 16 trapped ^'^^Yb"'" atoms. We 
directly control the amount of frustration by continuously tuning the range of interaction and directly measure 
spin correlation functions and their dynamics through spatially-resolved spin detection. We find a pronounced 
dependence of the magnetic order on the amount of frustration, and extract signatures of quantum coherence in 
the resulting phases. 



Cold atoms are ideal platforms for the simulation of frus- 
trated spin models, with the ability to tailor interactions 
with external fields and perform projective measurements of 
the individual spins [141- Neutral atomic systems are typi- 
cally limited to nearest neighbor interactions ifTTl . although 
geometrically-frustrated interactions can be realized in certain 
optical lattice geometries |20J. The natural long-range inter- 
action between cold atomic ions 11211 has led to the engineer- 
ing of Ising couplings between individual trapped ion spins 
ifTSI [191 122I I23I and the observation of spin frustration and 
quantum entanglement in the smallest system of three spins 
pTl. In this article we implement variable-range antiferro- 
magnetic (AFM) Ising interactions and transverse magnetic 
fields with up to iV = 16 atomic ion spins, using optical dipole 
forces. We directly measure the emergence and frustration 
of magnetic order through spatially-resolved imaging of the 
spins, in a system that approaches a complexity level where it 
becomes difficult or impossible to calculate the ground state 
order or the spin dynamics. 

We simulate the Ising model with long range antifeiTomag- 
netic interactions, given by the Hamiltonian {h = 1) 

H = Y,J.,ai^^<jy^-BY,al;\ (1) 

j<i i 

where a^^ (7 — x, y, z) are the spin-1/2 Pauli operators for 
the ith spin {i — 1,2,... N), B is the effective transverse 
magnetic field, and > is the Ising coupling matrix be- 
tween spins i and j, falling off with the lattice spacing |i — j| 



approximately as 

J^l ~ (2) 

where < a < 3 1.21 J . 

For B ^ Jij on all pairs, the spins are polarized along 
the effective transverse magnetic field in the ground state 
Itytj/ty • • • ) of '^he Hamiltonian in Eq. [T[ where Y\y) denotes 
a spin along the +y-direction of the Bloch sphere. As the ratio 
of B to the Ising couplings is reduced, the system crosses over 
to an ordered state dictated by the form of the Ising couplings, 
and the spectrum of energy levels depends on the range of the 
interactions. For any finite range interaction, the staggered 
AFM states |titi • ' ' ) and ||tlt • ' ' ) constitute the doubly 
degenerate ground state manifold at B = 0, with the degen- 
eracy arising from the time reversal or the global spin flip 
symmetry of the Hamiltonian. Here |t) and ||) are spins ori- 
ented along the Ising or x-direction of the Bloch sphere. Thus 
the system exhibits nearest-neighbor AFM or Neel ordering at 
sufficiently low temperatures. When the interactions are uni- 
form over all pairs of spins (a 0), the excitation gaps close 
significantly with maximum frustration, leading to a finite en- 
tropy density. In this case, any spin configuration with a net 
magnetization of zero (1/2) for even (odd) numbers of spins 
belongs to the ground state. 

Between the limits _B = and B ^ Jij, the energy spec- 
trum features a minimum gap, whose position and size de- 
pends on the amount of frustration in the system, or the inter- 
action range. (The interaction range is defined as the number 
of lattice spacings ^ where the interaction falls off to 20% of 
the nearest neighbor Ising coupling: ^ = 5^^".) Figure [la 
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Figure 1: Theoretical energy spectrum and critical gap in the long 
range antiferromagnetic Ising model (Eq. [TJ for A'^ = 10 spins, (a) 
For Ising exponent a — 1 (characteristic range of ^ = 5 sites), 
there are few low-lying energy states of the frustrated Hamiltonian 
as a function of the dimensionless parameter B/ Jo. The spacing be- 
tween the ground state and the first coupled excited state reaches a 
bottleneck at a critical value Be/ Jo with critical gap Ac. (b) Theo- 
retical dependence of Be/ Jo (dotted line) and Ac/ Jo (solid line) on 
the range of the interaction. As the interaction range increases, the 
competing long range couplings make it easier to create excitations 
and the critical gap is reduced, so a relatively small effective trans- 
verse field can break the spin ordering. Both parameters approach 
zero as a — >■ or ^ — >■ oo. Current experiments are performed with 
parameters in the shaded region. 



shows a few low lying energy states of the Hamiltonian in Eq. 
[T]for an interaction range of ^ = 5 (corresponding to a = 1). 
The first excited eigenstate merges with the ground state for 
small B/Jq and has the same spin order as the ground state 
near B/ Jq = 0. The critical gap Ac between the ground and 
the first coupled excited state determines the adiabaticity cri- 
terion ll24l . Figure [T|5 compares the position (red dotted line) 
and size (black solid line) of the critical gap of the Hamilto- 
nian for various ranges. As the range and hence the amount of 
frustration increases, the critical field is pushed towards zero, 
and the gap closes. The signature of frustration is the density 
of states near the ground state, and not merely the order in the 
ground state. In order to observe the effects of frustration, we 
therefore quench the system by ramping the effective trans- 
verse magnetic field faster than the critical gap (|_B/i3| > Ac) 
to populate the lowest coupled excited states. The observed 
spin order depends on the resulting degree of excitation and 
hence on the level of frustration. 

The spins are stored in a collection of ^^^Yb^ions con- 
fined in a linear radiofrequency (Paul) trap, with the effec- 
tive spin- 1/2 system represented by two hyperfine 'clock' 
states within each ion jfz) and Hz), separated by the hyper- 
fine frequency vhf = 12.642819 GHz |25|. The variable- 
range AFM Ising interactions are generated by applying off- 
resonant spin-dependent optical dipole forces |22| that drive 
stimulated Raman transitions between the spin states while 
modulating the Coulomb interaction between the ions in a 
controlled way (see Appendix I). The effective magnetic field 
is generated by simultaneously driving coherent transitions 
between the spin states with a 7r/2-phase shift with respect 
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Figure 2: Quantum phase transition from a paramagnet to an antifer- 
romagnet (Eq. [TJ, with Jij ~ |i — in a system of 10 spins, 

(a) Image of 10 trapped ions, with a distance of 22/im between the 
first and last ion. (b) Measured two-point correlation function be- 
tween a chosen spin (on the edge) and the other spins separated by 
r lattice sites, Ci,i+,. = (o-i^'cri^^''') — (cri^')(cri^^''), averaged 
over 4000 experiments for each value of the dimensionless parame- 
ter B/Jo- For 5/ Jo = 5 the spins are initially polarized along the 
transverse y-field with little correlation along the Ising x-direction. 
As the field is reduced the spins crossover to predominantly AFM 
states Ititi ■ ■ ■ ) ™d |itit ' ' resulting in alternating signs in the 
two point correlations with separation, (c) Measured occurance prob- 
ability of all 2^" — 1024 states at B/Jo ~ 5 (paramagnetic state, 
red trace) and 5/ Jo ~ 0.01 (AFM phase, black trace). The states 
are listed in binary order, with spin |4) = and |t) = 1- The resid- 
ual peaks in the red trace are consistent with detection errors biased 
toward states with many jt) spins such as 127, 255, and 1023. The 
two tall peaks in the black trace at 341 and 682 correspond to two 
Neel-ordered staggered AFM states, shown with camera images of 
these cases and contributing ~ 17% to the population. 



to the dipole force beams. At any time, we measure the state 
of the spins by illuminating the ions with resonant radiation 
and collecting state-dependent fluorescence on an imager with 
site-resolving optics [25 J . From this information we can ex- 
tract all spin correlation functions (see Appendix II). 

The quantum simulation begins by optically pumping all 
spins to the 14,2) state and then coherently rotating each spin 
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about the cc— axis of the Bloch sphere to initiaUze each spin 
in state |tj/) along the effective transverse magnetic field. The 
Hamiltonian (Eq. [T]) is then switched on with an initial field 
Bq « 5 Jo, where Jq is the average nearest neighbor Ising 
coupling, thus preparing the spins in the ground state of the 
initial Hamiltonian with a fidelity better than 97%. The effec- 
tive magnetic field is ramped down exponentially in time with 
time constant 400 /is (and no longer than 2.4 ms overall) to a 
final value B of the transverse field. At this point the Hamilto- 
nian is switched off, freezing the spins for measurement. We 
finally measure the x— or y— component of each spin (cri''') or 
(ay ^) by first rotating our measurement axes with an appro- 
priate global tt/2 pulse similar to the initialization procedure, 
before capturing the spin-dependent fluorescence on the im- 
ager. The experiments are repeated ~ 2000-4000 times to 
collect statistics on the resulting state. 

From these measurements, we can construct various order 
parameters appropriate for observing low energy AFM states. 
The staggered magnetization ~ X]i!=i(^l)*('''2:''') I de- 
lineates paramagnetic and AFM order, and also quantifies spin 
flip excitations. The fourth moment of this magnetization, 
known as a Binder cumulant = 3/2 — m^/2{m1)'^ varies 
from to 1 as the paramagnetic state gives way to AFM or- 
der and is scaled to remove finite size effects, and is aver- 
aged over independent experimental realizations. We also can 
form any correlation function of the spins such as the two- 
point correlation Ci_j — (ai'Vi/'') — (ctk "*) (cri"'''), allowing 
a direct probe of spin order for each experimental realiza- 
tion. The Fourier transform of this correlation function is 
the structure function S{k) = \ Y^^=i C{r)e^^'^\, where 

C{r) = I]m=i C'm,m+r IS flie average correlation over 
r sites in the chain. The structure function shows spin order 
versus wavenumber k, with S{k = it) singling out the pres- 
ence of the nearest neighbor Neel AFM order 

Figure |2] shows the onset of antiferromagnetic ordering 
in the quantum simulation of the frustrated transverse field 
Ising model in a system of 10 spins. Two-point spin cor- 
relations Ci r between a chosen edge spin and the other 
spins are presented in Fig. |2|3 at various stages in the ramp 
B/Jq — 5,0.25, and 0.01. For larger transverse magnetic 
fields, there are no appreciable correlations between the spin 
components along the Ising direction. As the ratio of B/ 
is lowered however, a zig-zag pattern emerges, with negative 
(positive) correlations between spins separated by odd (even) 
lattice spacings. For B/ Jq 0.01 the nearest-neighbor spin 
correlation reaches about 60%, and the correlation length (de- 
fined to be the distance at which the absolute correlation drops 
to 1/e of the nearest neighbor value) reaches about 6 lat- 
tice sites. The effective field was ramped exponentially down 
from B/ Jq « 5 with a time constant of 400 /its in this ex- 
periment. This ramping is not slow enough to be adiabatic, 
and the diabatic effects prevent the spin ordering to reach a 
perfect AFM phase. Figure [2]: shows the measured probabil- 
ities of all 2^" = 1024 possible spin states measured along 
the Ising x-direction, sorted in binary order with spin ||) = 
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Figure 3: Spin order vs the range of the Ising couplings for B/ Jo ~ 
0.01. (a) Structure function S(k) vs wavevector k for various ranges 
of AFM interactions. As the range of interaction increases, the grow- 
ing frustration suppresses the ground state order in the quantum sim- 
ulation, captured by the decreasing value of S{k — n). The er- 
ror bars include statistical fluctuations and estimated detection un- 
certainties. The detection errors are larger for the longest range of 
interactions, owing to spatial crosstalk from their closer spacing, (b) 
Distribution of observed states in the spin system, sorted according to 
their energy Ei that was previously calculated by diagonalizing Eq. 
[TJwith B = 0. Data is presented for two ranges (red for a — 1.05 
and blue for a = 0.76). The dashed lines indicate the cumulative 
energy distribution functions for these two ranges. 



and It) = 1. Since the detection fidelity of each spin within 
a chain is ^ 93%, the probability of correctly detecting a par- 
ticular state of the 10 spins is only 0.93^" « 48%. We recover 
effective detection efficiencies of ^ 98% per spin by post- 
filtering the measurements based on calibrating the known de- 
tection errors for each spin |26|, as described in Appendix 
II. The initial paramagnetic phase shown in red (B/Jq w 5) 
exhibits a roughly uniform probability of 1/1024 w 0.1% 
for each state (the residual peaks in the red trace are con- 
sistent with detection errors). The spin ordered phase shown 
in black (B/Jq ~ 0.01) displays the emergence of the two 
AFM states, each with an occupation probability of about 
8.5%. Other prominent peaks correspond to single spin-flips 
and other low-lying excitations from the two ground states. 

In Fig. [3]we probe the frustration in the system for various 
ranges of interactions. Here we look at the spin order achieved 
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Figure 4: Quantum coherence in tlie simulation. Probabilities of dif- 
ferent values of the total spin component along y in the initial po- 
larized state (red), when the transverse field is ramped to « O.OlJo 
(black), and when the transverse field is reversed back to its initial 
value (green). After reversal, the y-magnetization returns to ~ 63% 
of its initial value, indicating quantum coherence in the evolution. 
The trajectory of the transverse field (B, in green) and all the Ising 
couplings (J, in blue) is shown in the inset. 



in the quantum simulation when the external magnetic field is 
ramped down to B/ Jq w 0.01 for four different ranges of 
interactions. In Fig. [3^ we show the measured structure func- 
tion S{k) at wavevector values k ^ j^, j^,...,tt from the 
measured two-point correlation functions. To directly com- 
pare the different interaction ranges, we choose the same ex- 
ternal magnetic field ramp time constant, r = 0.4/Jo- As 
the range of interaction increases, the ground state AFM or- 
der (given by the structure function peak at A: = tt) disap- 
pears, denoting increased occupation of the excited states as 
the frustration grows. Figure ^ displays the observed dis- 
tribution of energy P{Ei) for a = 1.05 (shorter range) and 
a ~ 0.76 (longer range) power law exponents, along with 
the cumulative energy distribution function. The eigenener- 
gies of each configuration are calculated using Eq. [T] with 
B — {). Note that for the longer range interactions, excitations 
are more prevalent, and the energy gap between the ground 
and the first excited state is reduced, both signatures of in- 
creasing frustration in the system. The observed final entropy 
per particle S = —j^ J^i P{Ei)\ogP{Ei) is seen to increase 
as the interaction range grows, from S{a — 1.05) = 0.832 
to S{a = 0.67) = 0.903, which is also a signature of the 
increased frustration in the system. As a reference, the para- 
magnetic state distribution shows an entropy per particle of 
0.959, which is slightly less than unity due to detection errors. 

The above measurements of the state distribution concern 
only the diagonal components or populations of the density 
matrix. To characterize the quantum coherence in the simula- 
tion, we retrace the external magnetic field back to its initial 
value after ramping it down to almost zero, and measure each 
spin along the transverse {y) magnetic field. Figure |4] shows 
the distribution of the total transverse spin Sy = J^ii'^y^) 
at three different times: first at the beginning of the simula- 
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Figure 5: Spin order vs speed of ramp. The spins are initialized with 
B/ Jq = 5 and the transverse field is ramped exponentially down for 
6 time constants, and the experiment is repeated for various values of 
time constant r. Here we plot the scaled staggered Binder cumulant, 
Qs versus the total duration 6t of the ramp, and the solid line indi- 
cates the theoretical expectation of how the Binder order parameter 
should grow with time. As the Hamiltonian evolves more slowly, the 
observed spin order shows more ground state order, and less exci- 
tation for ramp times under ~ 2.5ms. However, for longer times, 
the spins become disordered, implying external decoherence in the 
system. 



tion, second when the transverse field has been ramped down 
to nearly zero, and third after the transverse field returns to 
its initial value. The initial state is ideally a delta function 
at Sy ~ 10, but finite detection efficiency broadens the dis- 
tribution. At the lowest value of the field (B/Jq sa 0.01), 
the transverse magnetization is distributed near Sy = 0, as 
the spins are presumably ordered along the Ising x-direction. 
When the external field is ramped back to its initial state, the 
distribution of the total spin returns toward the initial distri- 
bution, with a magnetization that is approximately 63% of its 
initial value, indicating that quantum coherence is maintained 
throughout the simulation. 

In order to probe decoherence in the simulation, we repeat 
the experiment with various ramping speeds of the effective 
magnetic field. In Fig. |5] we plot the AFM order parameter 
5s vs the total duration for the experiment for a long range 
coupling (a ~ 1.12 power law with distance) for iV = 10 
spins. Each data point represents the spin order achieved after 
ramping the magnetic field B down exponentially from 5 Jq 
for a total duration of 6 time constants. The AFM order grows 
with slower ramping, as expected, for up to r = 400/is. But 
we also observe a saturation and then decay in the spin order, 
which might indicate the presence of decoherence in the sys- 
tem at long times. During the simulation, spontaneous Raman 
scattering from the optical beams is expected to occur at a rate 
of less than lO^^fi ^ 6 s^^ per spin |,27J, which is consistent 
with separate measurements of the spin relaxation from a sin- 
gle spin and is therefore not expected to contribute to deco- 
herence given the timescales in the experiment. The phonon 
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Figure 6: Magnetic ordering in A'' = 16 spins, a. Image of 16 
trapped ions, with a distance of 30/xm between the first and last ion. 
b. Pair correlation function measured at various stages of the quan- 
tum simulation, for 5/ Jo — 5 (red) and B/ Jo — 0.01 (blue) in an 



AFM coupling falling off with distance as Jij 



among 



TV = 16 spins. While small amounts of staggered order are seen, it is 
tempered by the small gaps and frustration in the low energy states. 



c. In contrast, for all FM couplings (again with Jt- 



"^) the 



gaps are large and clear FM order is seen. Here the measured distribi- 
tion of magnetization is plotted. The paramagnetic phase of 16 spins 
is indicated in red, and after the field is ramped to nearly zero, the 
distribution clearly bifurcates, indicating population weighted heav- 
ily towards the FM states |m ■ • ■ ) and |ttt ■ ' ■ )• The resulting 
magnitude magnetization is approximately 73%. 



population is expected to be well under 10% for all the data 
presented here |28|. A principal source of decoherence ap- 
pears to be the intensity fluctuations in the Raman beams, due 
to beam pointing instabilities, and fluctuations in the optical 
power. 

Simulating adiabatic evolution in frustrated models is gen- 
erally slower than simulations without frustration, because the 
relevant energy gaps are smaller. As a direct comparison of 
these cases, we compare the quantum simulation of the long 
range AFM and ferromagnetic (FM) Ising models in a sys- 
tem of iV = 16 spins, shown in Fig. |6] For the FM exper- 
iment (Fig. [6};), we initialize the spins in the highest energy 
state with respect to the transverse field | iyiyiy • • • ) and 
ramp the field down as before. For the same simulation speed, 
we find that the best AFM nearest-neighbor correlation (Fig. 
^) is only ^ 30%, corresponding to a staggered magneti- 
zation of about 30%, while the simulation of the ferromag- 
netic model shows a clear ferromagnetic spin order across the 
chain, reaching ^ 73% magnetization. We have also observed 
a level of 74% FM magnetization emerging in = 18 spins. 
The cuiTent limit on calculating dynamics of fully connected 
spin models is approximately iV = 30 spins, a size that may 
be reached with technical upgrades in the hardware, including 
lower vacuum chamber pressures to prevent collisions with 
the background gas, better stability of the optical intensities, 
and higher optical power so that fluctuations in the beam in- 
homogeneities can be suppressed. 



Appendix I: Generating variable-range AFM Ising interactions 

The Ising interaction is generated by globally addressing 
the ions with two off-resonant laser beams at A = 355 nm 
EtII . intersecting at right angles with wavevector difference 
Afc along aprincipal axis of transverse ion motion |29|. These 
beams have beatnote frequencies vhf ^ IJ- which drive stim- 
ulated Raman transitions near the upper and lower motional 
sidebands of transverse motion in order to impart a spin- 
dependent optical dipole force li30ll . By setting the beatnotes 
sufficiently far from the sidebands, motional excitations can 
be made negligible, resulting in a nearly pure spin-spin cou- 
pling mediated by the Coulomb interaction ll22l . The effec- 
tive transverse magnetic field is generated by simultaneously 
driving a resonant stimulated Raman transition between the 
spin states with a beatnote frequency p and a phase that is 
shifted by 7r/2 with respect to the mean phase of the sideband 
fields. The resulting Ising coupling matrix Jij is given by a 
sum over contributions from each normal modes of collective 
motion at frequency Vm, 



Jij — ^"^VR / ^ 
m— 1 



(3) 



where i/r = h/M\^ = 18.5 kHz is the recoil frequency as- 
sociated with the dipole force, M is the mass of a single ion, 
&i „i is the orthonormal mode component of ion i with mode 
m, and Q, is the (uniform) single spin flip Rabi frequency, pro- 
portional to the laser intensity at each ion. The symmetric 
detuning /i of the beatnote from the spin-flip transition con- 
trols the sign and range of the interactions [22 J. When /i is set 
larger than the highest (center-of-mass or COM) mode fre- 
quency vi, every interaction is AFM, and we can empirically 
approximate Eq. |3]as falling off with distance as a power law 
J,j w Jo/|i-j|" withO < a < 3 |21| and Jo « 1/iV. While 
the COM mode mediates a uniform interaction between all 
pairs of spins, the other modes introduce non-uniformity in the 
interactions, and effectively reduce the range of AFM interac- 
tion. In practice, we control the interaction range by changing 
the bandwidth of the transverse mode spectiTim, achieved by 
varying the axial confinement of the ions in the Paul trap. The 
Ising couplings J^ depend not only on the spatial separation 
\i — i\, but also on the site i itself due to the finite size of the 
system, with ~ 10% inhomogeneities across the chain. We 
average over all the couplings between spins separated by a 
given number of lattice sites to estimate the power law range 
exponent a in Eq. |2] 1 19J . 

In the experiment we use global Raman beams each with 
an optical power of ^ IW, having horizontal and vertical 
waists of ^ 150 /^m and ~ 7 /im respectively to address 
the ions. This produces a spin-flip Rabi frequency Vl ^ 600 
kHz on resonance, with less than 5% inhomogeneity across 
the chain. We set the beatnote detuning to ^ « i^i + 377^!, 
where r/ — yjvji/vi is the single ion Lamb-Dicke parameter. 
This keeps the (primarily COM) phonon excitation probabil- 
ities sufficiently low for any setting of the range. The typical 
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nearest neighbor Ising coupling is Jq ~ 1 kHz for iV = 10 
spins. In principle, the Ising interaction range can be varied 
from uniform to dipolar (0 < a < 3), but in this experi- 
ment the axial frequencies was only varied between 0.62 MHz 
and 0.95 MHz, and given the COM transverse frequency of 
vi =4.1 MHz, this results in a range of Ising power-law ex- 
ponents 0.7 < a < 1.2, or a variation of the range of interac- 
tions between ^ = 4to^ = 10 sites. 
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Appendix II: Accounting for finite detection efficiency of 
N-particle correlations 

We detect the spin states using spin dependent fluores- 
cence collected through f/2. 1 optics on an intensified charge- 
coupled-device (ICCD) camera or a photomultiplier tube 
(PMT). The spin state \ fz) fluoresces from the near resonant 
detection beam, and appears bright, while the spin state 
scatters little from the off-resonant detection beam, thus ap- 
pearing dark. The imager has single site resolution, allow- 
ing us to directly measure the two point correlations to probe 
the AFM order The single spin detection efficiency is about 
e = 93% on the ICCD imager, and about 98% on the PMT, 
the reduced efficiency on the ICCD being due to electronic 
and readout noise. To account for the spatial overlap of the 
fluorescence from neighboring bright ions, we fit each single 
shot image in the experiment to a sum of N Gaussians, where 
N is the number of ions. The center and width of the Gaus- 
sians are pre-calibrated from images of all spins prepared in 
the bright states, with background subtraction from all spins 
prepared in the dark states. The probability of correctly iden- 
tifying a = 10-body spin state is only sa 48%, so we 
post process the detected states to account for this finite de- 
tection efficiency |26|. The probability of incorrectly assign- 
ing an A^-qubit state \i) to the actual underlying state |j) is 
Mij = (1 — g)^>je^^ftj , where f3ij is the number of posi- 
tions that the A^-qubit state \j) differs from \i) through bit 
flips. The observed probability distribution of all 2^ states 
is given by P/ = J^j ^^ijPj^ where Pj is the underlying ac- 
tual distribution of states, which can be obtained by simply 
inverting the matrix Mij and forming Pj ~ M^^Pj. This 
increases the effective detection fidelity to about 98%, equiv- 
alent to that with the PMT. Some entries of the post-processed 
probabilities are slightly negative, due to fluctuations in ab- 
solute fluorescence levels that impact the values in the matrix 
Mij during measurement. 
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